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ABSTRACT 

We present a spherically symmetric model for the origin and evolution of the temperature 
profiles in the hot plasma filling galaxy groups and clusters. We find that the gas in 
clusters is generically not isothermal, and that the temperature declines with radius 
at large distances from the cluster center (outside the core- and scale radii). This 
temperature profile is determined by the accretion history of the halo, and is not 
quantitatively well-described by a polytropic model. We explain quantitatively how the 
large-scale temperature gradient persists in spite of thermal conduction and convection. 
These results are a consequence of the cosmological assembly of clusters and cannot 
be reproduced with non-cosmological simulations of isolated halos. We show that the 
variation in halo assembly histories produces a ~ 10% scatter in temperature at fixed 
mass. On top of this scatter, conduction decreases the temperature of the gas near the 
scale radius in massive clusters, which may bias hydrostatic mass estimates inferred 
from x-ray and SZ observations. As an example application of our model profiles, we 
use mixing-length theory to estimate the turbulent pressure support created by the 
magnetothermal instability (mti): in agreement with our earlier mhd simulations, we find 
that the convection produced by the mti can provide 5% non-thermal pressure support 
near rsoo- The magnitude of this turbulent pressure support is likely to be non-monotonic 
in halo mass, peaking in ^ lO^"*'^ Mq halos. 
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1 INTRODUCTION 

X-ray observations of the hot, diffuse gas in galaxy groups and 
clusters suggest negative temperature gradients at large radii 
(Leccardi & Molendi 2008; George et al. 2009; Simionescu 
et al. 2011). This observation is somewhat surprising because 
negative temperature gradients are susceptible to a convec- 
tive instability known as the magnetothermal instability, or 
MTI (Balbus 2001). Without a clear source of free energy to 
maintain the convection, it seems unusual that clusters should 
so uniformly be found in unstable states. This observation 
becomes even more surprising in clusters more massive than 
~1O^^'^M0, where the timescale for heat to diffuse through 
the ICM can be much shorter than the Hubble time. Left to 
their own devices, both conduction and convection tend to 
erase temperature gradients, and one might expect them to 
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make the gas isothermal. In this paper, we study how the 
assembly of clusters creates large-scale temperature gradi- 
ents and maintains them in spite of convection and thermal 
conduction. 

Our subject is not purely academic. Systematic trends 
in temperature gradients with mass and redshift may influ- 
ence the conversion of observable quantities (such as the x-ray 
surface brightness) to thermodynamic quantities (gas density, 
pressure, etc.). Understanding the origin of the temperature 
gradients in clusters would also enable us to calculate how the 
turbulence produced by the mti depends on halo mass or red- 
shift. These trends create systematic variations in non-thermal 
pressure support and might affect current efforts to use the 
cluster mass function to constrain cosmology (see Parrish et al. 
2012, and references therein). Thus, while the temperature 
profiles in galaxy clusters present an interesting puzzle in their 
own right, understanding the processes which control them 
may also find useful application in cluster cosmology. 
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Several studies, including Dolag et al. (2004), Burns et al. 
(2010), and Ruszkowski et al. (2011) have begun to address the 
effects of conduction and the mti in the ICM using cosmologi- 
cal simulations. These results are computationally expensive, 
however, and can be difficult to interpret. For example, the 
MTI is expected to be a sub-dominant, but significant, source 
of turbulence in clusters. Using cosmological simulations to 
study its trends with mass and redshift would require very 
careful calibration of other sources of turbulence. Thus, we 
feel that a simplified treatment which affords an intuitive 
understanding of the results remains useful. 

The temperature and entropy profiles in the ICM are 
related by hydrostatic equilibrium. Accordingly, this paper 
closely follows earlier work by Tozzi & Norman (2001) and 
Voit et al. (2003), who study the entropy profiles in clusters. 
However, the processes controlling temperature gradients in 
clusters are slightly more subtle than those which determine 
their entropy profiles. Radial variations in entropy tend to be 
much larger than those in temperature, so small differences 
in the entropy profiles translate to much larger differences in 
temperature profiles. For the same reason, thermal conduc- 
tion has a more pronounced effect on the temperature profile 
than on entropy and we must include it in our analysis. Thus, 
despite significant similarities to both Tozzi & Norman (2001) 
and Voit et al. (2003), our models represent a generalization of 
these earlier studies and we use them to explore different astro- 
physical applications. Komatsu & Seljak (2001) also present 
analytic models for temperature profiles in clusters. As we 
describe below, however, our method differs significantly from 
theirs because we do not assume a polytropic model for the 
gas; instead we directly calculate its thermodynamic state. 

We describe our general methodology in the next section, 
followed by two sets of models. We begin in section 3 with 
an idealized but intuitive model. This section illustrates the 
key process that sets temperature gradients in clusters, but is 
too idealized to be directly compared with real clusters. We 
generalize this model in section 4 to more accurately reflect the 
formation histories and gravitational potentials of clusters. We 
also study the influence of thermal conduction on our results. 
We close in § 5 with a brief summary of our conclusions and a 
description of our future plans to apply them to observations 
and to more realistic simulations. 



2 METHOD 

Following Tozzi & Norman (2001) and Voit et al. (2003) (here- 
after V03), we model the cosmological assembly of a galaxy 
cluster and use the evolving properties of its accretion shock 
to calculate the thermal state of the gas in its ICM. This differs 
from the approach taken by Komatsu & Seljak (2001), who 
adopt a polytropic equation of state and thus assume a ther- 
modynamic state for the gas. Since the dynamical timescale in 
clusters is typically shorter than the timescales for either ac- 
cretion or thermal conduction, we approximate the dynamics 
of the gas by assuming that it is in hydrostatic equilibrium at 
all times. We focus our attention on the temperature profile 
generated by the accretion history of the halo and on how 
this profile is modified by thermal conduction. As discussed 
in section 1, we build on previous work that has focused on 
the origin of cluster entropy profiles. 

Galaxy cluster halos are tri-axial and grow partly by ac- 



creting smaller sub-halos. The extent to which these properties 
influence the temperature profiles in clusters is not entirely 
clear; including them in our analysis would require cosmo- 
logical simulations, however, and would likely obscure our 
results. Instead, we model clusters as spherically symmetric 
and we assume that they grow primarily via smooth accre- 
tion. Though these approximations are not entirely realistic, 
they are partially justified in V03, who compare models based 
on smooth accretion with more realistic ones based on hier- 
archical structure formation: the differences are modest for 
high-mass halos and mostly reflect changes in normalization 
due to clumping factors. 

Each shell of material accreted by the cluster contains 
both dark and baryonic matter. We do not solve for the 
evolution of the dark matter; instead, we assume that the dark 
matter virializes rapidly and we parametrize it with a time- 
dependent fit to the potential (we neglect the gravity due to the 
baryons). We consider both simplified, isothermal potentials 
(§ 3) and more realistic fits to cosmological simulations (§ 4). 

The baryons enter the cluster through a shock with a 
density pb,i and velocity Vi determined by the accretion rate 
and by the depth of the potential: 

/b 1 9M,h 



pb,i Vi 



4^ '^'h 



dt 



In the above, r^i-i is the radius of the accretion shock, Mgh is 
the total mass contained within it, and /b ~ 0.17 is the cosmic 
baryon fraction. The parameter ^ = rgh/^ta parametrizes the 
strength of the virial shock, and rta is the "turnaround" radius, 
at which the kinetic energy of the shell vanishes (Gunn & Gott 
1972).^ We assume that the IGM is cold enough that the gas 
enters the cluster through a strong shock (with an upstream 
Mach number M ^ 1). Thus, the post-shock density and 
pressure are given by: 



fsh = |Pb,i 1 T^Vi 
Psh = 4/9b,i, 

where the factor of 4/3 in front of the velocity transforms the 
infall velocity from the frame of the cluster to the frame of 
the shock. Thus, the post-shock entropy (K = k-BTn~^^^) is: 



-2/3 



( a In M.h \ ^ 
. \ 9' ) 



(1-0' 



1/3 



(1) 



Equation 1, above, is identical to equation 8 in V03. As empha- 
sized in V03, apart from an overall normalization oc {M / f'aY^^ , 
the entropy profile depends only on the accretion history and 
on the shock strength ^. 

The jump conditions above represent the conservation 
of mass, momentum, and energy and thus refiect the total 
pressure behind the shock. In general, however, the electron 
and ion pressures may differ. The difference may be significant 

^ Note that identifying 5 with r^jj/rta assumes that the baryons 
do not separate from dark matter before reaching the virial shock. 
If this were not the case, the baryons would feel a time-dependent 
potential due to the dark matter. Tozzi & Norman (2001) find that 
this may introduce a ~10% correction to the infall velocity. 
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in clusters because the shock preferentiaUy heats ions and be- 
cause the timescale for electrons and ions to equilibrate is long 
(Fox & Loeb 1997). We don't distinguish between electron and 
ion temperatures in our calculation because the temperature 
difference does not effect hydrostatic equilibrium and thus 
should not influence our solutions. Moreover, the simulations 
by Rudd & Nagai (2009) show that this temperature difference 
is modest within the virial radius of the cluster. 

Our assumption of a strong virial shock may not be valid 
for the innermost shells of material, which accreted when the 
IGM was hotter, and when the gravitational potential of the 
halo was shallower, than they are today. Hence, our model 
will not accurately reproduce the gas profiles near the centers 
of clusters (Tozzi & Norman 2001). Many other processes, 
including cooling, heating, and the formation of the central 
galaxy also affect the structure of the core, however (Voit 
2011). We instead focus on the gas at larger radii, near the 
virial radius. Both the large mass and the long cooling time 
of this gas enable us to ignore galaxy formation processes at 
smaller radii (e.g. Voit & Ponman 2003; Voit 2011). 

After undergoing the virial shock, the gas entropy evolves 
due to thermal conduction. Thermal conduction is highly 
anisotropic in the ICM because the electron mean free path 
is much longer than its gyroradius. Anisotropic conduction 
renders the ICM unstable to the mti (Balbus 2001), how- 
ever, which generates turbulence and may isotropize the mag- 
netic field (McCourt et al. 2011; Parrish et al. 2012). We 
therefore parametrize thermal conduction through an effective 
isotropic conductivity. We introduce the effective conductivity 
Kofi = /spKo, where /sp ~ 1/3 is a suppression factor due 
to the magnetic field. ^ Thus, the evolution of the entropy is 
determined by the following energy equation: 



dt 



3P 



vr 



(2) 



In the simplifying case that ^cff — > 0, the entropy of each shell 
is a constant in time and equal to -K'ah, given by equation 1. 
In addition to the conductivity k, we also use the thermal 
diffusion coefficient Xc = i^a/ineks), which has units of cm^/s. 

As mentioned above, anisotropic thermal conduction in 
the ICM drives a convective instability known as the MTI. This 
convection carries a heat flux which should technically be 
added to equation 2. However, the convective heat flux is 
smaller than the conductive flux by a factor of 



where ttIo and mi are the electron and ion masses, H is the 
pressure scale-height, Ao is the electron mean-free-path, and a 
is the mixing-length parameter of the convection. This ratio is 
small (~ 10"^) in the ICM, enabling us to ignore the convective 
heat flux in equation 2 (cf. Parrish et al. 2008). 



^ Though this approximation is suitable for our purposes, we note 
in passing that the use of an isotropic conductivity significantly 
alters other processes in the ICM, such as convection (Balbus 2001; 
Quataert 2008) and thermal instability (Sharma et al. 2010; Mc- 
Court et al. 2012; Sharma et al. 2012), and is thus not suitable for 
more detailed dynamical studies. Interestingly, a suppression factor 
/Sp ~ 1/3 turns out to be appropriate even if the magnetic field is 
tangled on very small scales (Narayan & Medvedev 2001). 



Our spherically symmetric model is most easily repre- 
sented in Lagrangian coordinates with the enclosed mass as 
the independent variable. Thus, the continuity equation and 
the equation for hydrostatic equilibrium take the form: 



dr 
dM^ 
dP 



1 



47J-J-2 

9 



3/5 



47rr2 



(3) 
(4) 



where Mb is the baryonic mass contained within the radius 
r, and we adopt the entropy K and pressure P as our pri- 
mary thermodynamic variables. In this initial study, we ignore 
sources of non-thermal pressure support (such as cosmic rays, 
magnetic fields, or turbulence). However, we discuss in sec- 
tion 5 our plans to self-consistently include turbulence gener- 
ated by the mti and to compare these results with cosmological 
simulations. 

Equations 1-4 completely specify our model except for 
the mass accretion history Msh{t) and the gravitational field 
g{r, t) of the halo. We assume that these are determined by 
cosmology and are unaffected by the baryonic formation of 
the cluster. We describe a simple, idealized model for the halo 
formation in section 3 and a more realistic model in section 4. 



3 SIMPLIFIED ADIABATIC MODELS 

Before studying the consequences of thermal conduction, it is 
useful to understand the 'baseline' temperature profile gener- 
ated by the halo's mass accretion history. Thus, we begin with 
adiabatic models which ignore thermal conduction. We further 
isolate the influence of the accretion history by assuming the 
gravitational potential is isothermal and that the accretion 
rate of the cluster is independent of time: Msh(i) = Mq x t/to. 
As we show later, this prescription contains only a single free 
parameter and is perhaps the simplest nontrivial model of the 
process we wish to study. Though the results in this section 
cannot be directly applied to clusters, they highlight some of 
the key physics determining the temperature gradients in clus- 
ters, and will assist in our interpretation of the more detailed 
models in section 4. 

The next step in our model is to determine the shock 
strength ^. The most logical choice would be to calculate 
^(Mb) so that the cluster was in hydrostatic equilibrium at 
every epoch; in fact, this is formally required to use equation 1 
for the post-shock entropy. In this section, however, we make 
the simplifying assumption that ^(Afb) is a constant. Thus, 
the models presented in this section are not entirely self- 
consistent. Our goal in this section is only to obtain an intuitive 
understanding of how the accretion rate of a halo influences 
its temperature gradient. We present more detailed models, 
with more accurate results, in section 4. 

Since we have assumed a solution for ^(A/b) and that the 
evolution of the gas is adiabatic, the entropy profile J('(Afb) 
is uniquely determined at all times by equation 1. Thus, it 
suffices to solve hydrostatic equilibrium only at the present 
epoch - this solution cannot depend on the state of the cluster 
at earlier times. Hence, we need not track the evolution of the 
cluster, and the equations determining the state of the gas 
reduce to ordinary differential equations. 
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3.1 Method 

Before solving the equations of our model, we recast them in a 
more versatile dimensionless form. Since we have a system of 
ordinary differential equations with the shock as one boundary, 
we de-dimensionalize the equations in this section using the 
properties at the shock radius. (Note that this differs from the 
usual convention of using the virial radius of the underlying 
dark- matter potential.) We introduce the constants i?sh and 
Mo, which represent the shock radius of the cluster at the 
present time to and the total mass enclosed within it. (In 
an isothermal potential, Rsh and Mo each differ from the 
virial radius and mass by a factor of 2^.) We also define a 
dynamical time tdyn = (GMo/R^-^^)^^^^ , along with the spatial 
coordinate x = r/Rsh and the Lagrangian mass coordinate 
j/ = Mb//bMo. 

We introduce the dimensionless gas variables Pi, pi, and 
Ki via: 

— ^ X pi 



Pb 



4- RL 
/b gmS 

47r Rf, 



X Pi 



h 
5/3 



(5a) 

(5b) 
(5c) 



Thus, the equations for mass conservation and hydrostatic 
equilibrium become: 



dy 33/5 1^ to 

dPi _ _^ 

dy x^ ' 



-2/5 4/S 



y 



x'^Pl 



3/5 



(6a) 



(6b) 



which we solve subject to the boundary conditions at the 
shock: 



x{y = 1) = 1 

A(j/ = l) = ^(l-C)^/^(^ 



and at the center: 



c(j/ = 0) = 0. 



(7a) 
(7b) 

(7c) 



With these definitions, the solution is independent of the 
parameters /b, p, Rsh, and Mo. Furthermore, the system 
is over-determined with three boundary conditions and two 
equations. The shock radius ^ is therefore an eigenvalue which 
must be chosen to meet the inner boundary condition in 
equation 7c. 

Thus, the only free parameter in our system of equations 
is the ratio tdyn/to, and we expect to find a one-dimensional 
family of models. Since the average density of any dark matter 
halo p ~ 200pcrit is independent of mass in hierarchical struc- 
ture formation, we do not expect the dynamical timescale tdyn 
to vary strongly among clusters at any cosmological epoch. 
The ratio tdyn/to thus measures the age (or, equivalently, the 
assembly rate) of the halo. 

We note that the following solution to equations 6 and 7: 



Pi = 



2a;2 



x = y, C = 



(8) 



exists when the assembly rate satisfies: 

^dyn 



We describe the physical significance of this special, isothermal 
solution in section 3.2. 

For other values of the assembly rate tdyn/ioi we solve 
equations 6 and 7 numerically using a predictor-corrector 
method, and we solve the eigenvalue problem for ^ with a 
bisection search. We avoid the singularity in equations 6a 
and 6b by solving them on a logarithmic grid. Since we cannot 
apply the boundary condition equation 7c in the logarithmic 
coordinates, we obtain an approximate boundary condition at 
a finite radius by expanding the equation for hydrostatic equi- 
librium near a; = 0. Assuming that the temperature remains 
finite, this equation becomes: 



dlnp 



-1/2 



ci In r 



Combined with the fact that K oc Tp' 
that p (X r~^, and thus that 

T Tvir, and 



-2/3 



oc y 



y 



^Z'', this implies 

(9) 
(10) 



as X 0. We use equation 10 as a boundary condition for 
the eigenvalue problem and equation 9 to check the accuracy 
of our integration. The analytic solution in equation 8 also 
permits a more pedantic test of our method. 



3.2 Results 

Figure 1 shows the solution for the dimensionless shock radius 
£,{tdyn/to) from our numerical calculations. In the limit of 
very slow accretion (i.e. as tdyn/to — 0), the shock is com- 
paratively weak and the dimensionless shock radius ^ — >■ 1. 
As the assembly rate tdyn/to increases, the shock radius ^ 
decreases monotonically. Thus, the virial shock moves inwards 
as the accretion rate increases. This result seems intuitively 
reasonable, as a higher accretion rate implies a higher ram 
pressure behind the shock. A stronger shock (or smaller 5) is 
thus required to hold back the infalling material and to keep 
the ICM in hydrostatic equilibrium. 

Interestingly, figure 1 indicates a maximum assembly 
rate around tdyn/to ~ 0.39 at which ^ — >■ 0. Beyond this 
point, thermal pressure alone cannot hold back the accretion 
shock and hydrostatic equilibrium becomes impossible. This 
represents an extremely rapid accretion rate, however, with the 
halo forming over only ~ 2.5 dynamical times; our quasi-static 
model for the ICM breaks down in this limit. A fully dynamical 
calculation (e. g. a simulation) with such a high accretion 
rate would likely produce a model ICM with significant time 
dependence and turbulent pressure support, but with a finite 
shock radius. 

Figure 1 shows the isothermal solution (eq. 8) as a blue 
square and a point representative of a massive cluster [with 
tdyn = O.lHo^ and to ~ {2/3)Ho^] as a red triangle. Our 
model predicts that ^ ~ 0.6 for this fiducial cluster, similar 
to what has been expected in the past (e. g. Rees & Ostriker 
1977). Note that the isothermal solution requires a much faster 
assembly than is typical for galaxy clusters; this simple model 
thus suggests that clusters should not be isothermal. 

Figure 1 directly illustrates the effect of the accretion 
rate on the location of the virial shock. The gas properties 
in the ICM must match onto jump conditions at the shock; 
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o 




tdyn/to 

Figure 1 . Shock radius ^ as a function of the assembly parame- 
ter tdyn/^O for the simphfied models described in section 3. The 
timescale to = M/M is the age of the halo. Points further to the 
left on this plot correspond to clusters which form slowly and points 
further to the right correspond to clusters which form rapidly. The 
red triangle represents the approximate value of tdyn/io expected 
for cluster halos. The blue square marks the isothermal solution in 
equation 8. Clusters in the white region of the plot have negative 
temperature gradients, clusters in the pink region have positive 
temperature gradients, and clusters in the gray region cannot ex- 
ist in our steady-state model. Real clusters are expected to have 
moderate, negative temperature gradients based on this analysis 
(see figure 3). The black line illustrates the calculation with an 
isothermal potential, and the gray line shows results with a more 
realistic NFW potential. 

thus, by moving the shock radius, the accretion rate also 
influences the thermodynamic structure of the ICM. Figure 2 
demonstrates this by showing profiles of the gas density and 
entropy for models with different assembly rates tdyn/^o- The 
profiles are nearly isothermal, with p ~ and K ~ r^^'', 
and are broadly consistent with determinations from x-ray 
data (Croston et al. 2008; Cavagnolo et al. 2009). 

The small deviations from power-laws in the density and 
entropy profiles lead to significant temperature gradients, how- 
ever. We show this explicitly in figure 3, where we plot the 
temperature profiles for the models from figure 2. The temper- 
ature profiles are nearly linear, and gradients of either sign are 
possible, depending on the assembly rate tdyn/io- This result 
is not surprising, since the central temperature must equal 
the virial temperature of the halo (eq. 9) and the temperature 
at the shock is dictated by the jump conditions. A roughly 
linear interpolation between these boundary conditions seems 
reasonable given the simplicity of the model. The isothermal 
solution with tdyn/^o = \/3/32 divides models with negative 
and positive temperature gradients. In reality, most clusters 
satisfy tdyn/to ~ 0.1 — 0.2; in this case our model predicts 
temperature profiles which decrease by a factor of ~ 2 from 
the center to the shock radius. This result is in line with re- 
cent observations of the gas temperature near the virial radius 
(George et al. 2009; Simionescu et al. 2011). 



I I I I I I I r 

trlyn/U = 0.01 




I I I I I I I l_ 

10-1 1 



Figure 2. Profiles of density and entropy for representative mod- 
els from fig. 1 with different values of the assembly parameter 
tdyn/to. Thick gray lines show typical power-law slopes derived 
from x-ray observations (Croston et al. 2008; Cavagnolo et al. 2009) 
(the normalization is arbitrary). Massive clusters should mostly lie 
between the solid blue and long-dashed red lines. Our solutions are 
approximately, but not exactly, power-laws. The deviations from 
power-law behavior are dictated by the outer boundary condition 
on the pressure and determine the temperature profile (fig. 3). This 
boundary condition depends on the ram pressure behind the shock, 
and thus on the speed with which the cluster formed (parametrized 
by tdyn/to in this model). As in fig. 1, thick curves show models 
with simplified, isothermal potentials and thin curves show more 
realistic models with NFW potentials. 



Thus, while we solve the system of equations 6a and 6b 
and boundary conditions 7a-7c for the shock radius and for 
the structure of the ICM, we find that the boundary conditions 
essentially dictate the temperature profile. The temperature 
must reach the virial temperature of the potential at the center 
of the halo, and it must match onto the jump conditions at the 
shock. The solution to hydrostatic equilibrium then implies 
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idyn/tn = 0.2 

tdyn/to = 0.25 

tdy„/t„ = ^3/32 
tiyjta = 0.37.5 



0.0 



0.2 



0.4 0.6 



0.8 



1.0 



Figure 3. Temperature profiles for representative models from 
fig. 1 with different values of the assembly parameter tdyn/to- The 
black line corresponds to the isothermal solution (eq. 8); massive 
clusters should lie between the blue and red lines. Our model implies 
temperature profiles which decrease by a factor of ^ 2 from the 
center to the virial radius, in line with x-ray observations. Thin, 
light lines show the results of calculations with NFW potentials, 
rather than isothermal ones (see § i for details). The temperature 
gradients are similar outside the scale radius of the halo. 



a nearly linear interpolation between these two boundary 
conditions. 

As presented here, the temperature profiles in figure 3 
may seem specific to our assumption of an isothermal potential. 
To demonstrate that this is not the case, we have included 
calculations with nfw potentials in figure 3 (shown as thin, 
light lines). Although the dilTerent potential has a dramatic 
effect on the temperature profile within the scale radius, the 
trend between the assembly rate and the overall temperature 
gradient at large radii is similar. We include these lines only 
for illustration, but present much more detailed models with 
NFW potentials in the following section 4. 

The simple model presented in this section suggests that 
the assembly rate of the halo (or, equivalently, the ram pressure 
behind the accretion shock) dictates the large-scale temper- 
ature gradient in the ICM. The assembly rates of massive 
clusters are such that they should have moderate, negative 
temperature gradients outside the scale radius. This is one of 
our primary findings. In what follows, we show that this result 
holds true even when we relax the simplifying assumptions in 
this section. We also study how thermal conduction modifies 
this 'baseline' temperature profile. 



4 CONDUCTION AND REALISTIC ASSEMBLY 
HISTORIES 

Our models with isothermal potentials and linear accretion 
histories are especially transparent. The models are very ide- 
alized, however, and it is not clear how accurately they carry 



over to real clusters. In this section, we generalize our results 
to include more realistic potentials and accretion histories; we 
also include thermal conduction in our calculation and study 
its influence on the temperature profiles in clusters. 



4.1 Method 

4.1.1 Coordinates and AssumpUons 

Introducing realistic potentials and accretion histories into 
our model necessitates a few changes to our method. Though 
the unit system introduced in section 3 is ideal for our model 
equations, the virial mass of the halo (as opposed to the mass 
enclosed by the virial shock) is an eigenvalue of the problem 
and cannot be specified ahead of time. Since we want to study 
the variation in temperature profiles at fixed virial mass, we 
must alter the equations and boundary conditions slightly. 

We adopt a unit system based on the virial mass Mvir, 
and radius rvir,o of the halo at redshift zero. Following our 
approach in the last section, we introduce the dimensionless 
variables Pi, pi, and Ki via: 



Ph 



rl-o 



X pi 



47r r^ 

Pi 
5/3 ■ 
Pi 



(11a) 
(lib) 
(11c) 



We also define the spatial coordinate x = r/r^ii, the La- 
grangian coordinate y = Mb//bAfvir, and the dynamical time 
tdyn = (GMvir/^vir)"^^^- Note that our definitions of x, y, 
and tdyn use instantaneous values of r^ir and A/vir, while our 
definitions of pi, Pi, and Ki are normalized to rvir, and 
A/vir, 0- Thermodynamic quantities in our calculation (e.g. 
Ki) are thus directly comparable at different redshifts, while 
coordinates (e. g. x) are not. 

We define m = Mvir{t)/Mvir,o, which tracks the forma- 
tion of the halo and functions as a time coordinate. We take the 
virial radius to be r2oo, the radius within which the mean den- 
sity of the halo is 200 times the critical density of the universe. 
Thus, tdyn ~ {10H)~^ , where H is the Hubble parameter. 



4.1.2 Dark Matter 

As discussed in section 2, we do not solve for the evolution 
of the dark matter; instead, we assume that the dark matter 
evolves independently of the baryons and we parametrize 
it using fits to halos from cosmological n-body simulations. 
In particular, we assume that the dark matter follows an 
NFW distribution at all times (Navarro et al. 1997), with a 
constant concentration parameter c = 5, as is appropriate for 
actively forming, massive halos (Zhao et al. 2009). The only 
free parameter in this model for the dark matter is the mass 
accretion history of the halo. 

We use fits to halo mass accretion histories of the form 



m{z) = [(1 + z)*" exp(- 



(12) 



derived by McBride et al. (2009) from the Millennium simula- 
tion (Springel et al. 2005). This fit is calibrated to the "friends 
of friends" mass AfpoF (Davis et al. 1985), which is similar to 
our choice of M200 for the virial mass (White 2001). Figure 4 
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Figure 4. Parameter space for dark matter accretion histories 
based on the Millennium Simulation. The exponents 7 and b are 
defined in equation 12. Colors demarcate the different accretion 
"types" from McBride et al. (2009) and gray lines show contours of 
the formation redshift 2f such that Mvir(zf) = 0.5Mvii{z = 0). At 
redshift zero, Type III halos are accreting very slowly and Type 
halos are accreting rapidly. The red curves mark contours of the 
probability density for IO^^Mq clusters to have a given accretion 
history (taken from the Appendix of McBride et al. 2009). The 
contours for IO^'^Mq halos are very similar to those shown here. 



shows the parameter space for the exponents 7 and b, along 
with the definitions of Type I, II, III, and IV accretion histories 
from McBride et al. (2009).^ For reference, the red curves in 
figure 4 show contours of the probability density function for 
halo accretion parameters from the appendix of McBride et al. 
(2009). We also show contours of the "formation redshift" Zf 
at which the cluster reaches half its present mass; most of the 
halos formed at a redshift between 0.5 and 1. 

Note that a Type I accretion history is roughly exponen- 
tial in redshift, which is typical in ACDM (Wechsler et al. 
2002). At fixed mass. Type halos are younger, while Type II 
and Type III halos are older, than Type I halos. A Type IV 
accretion history corresponds to mass loss at late times; this 
would imply a negative ram pressure at the virial shock in our 
quasi-equilibrium model. These cases should be studied with 
fully dynamical, cosmological simulations. Fortunately, at the 
high masses we wish to study, a relatively small fraction of 
the total halos exhibit Type IV accretion histories. Moreover, 
these are systems which have recently undergone major merg- 
ers; they are likely to be morphologically disturbed and may 
be excluded from cosmological samples. 

In order to calculate the strength of the virial shock, 
we require an estimate for the turnaround radius rta (§ 2). 
Unfortunately, fits to rta from n-body simulations do not seem 

^ McBride et al. (2009) do not distinguish between Type II halos 
with negative and positive values of the exponent b. Since this 
distinction is important in our application, we denote Type II halos 
with b < as Type 0. 



to be available. Therefore, we simply use the virial theorem 
to estimate that the turnaround radius is twice the virial 
radius. The shock radius is then given in terms of the virial 
radius by Vsh ~ 2^rvii. We note that this approximation 
may over-estimate the turn-around radius (cf. Diemand et al. 
2007), causing us to predict temperature profiles which are 
too shallow. This is one of the primary sources of uncertainty 
in our models. 

The quasi-equilibrium model for the dark matter de- 
scribed in this section greatly simplifies our method by elimi- 
nating the need to solve for the dark matter dynamics (e. g. by 
using an n-body simulation or by solving the Jeans equations) . 
This model is is ambiguous outside the virial radius, however. 
Since the virial shock typically lies exterior to the virial radius, 
it is not a priori clear what value of M^h to use in equation 1 
for the post-shock entropy. We proceed by presuming that the 
gas and dark matter first separate at the virial shock; thus, 
the gravitating mass Msh in equation 1 corresponds to the 
mass Mvir. After the shock, the gas remains in hydrostatic 
equilibrium at rsh while the dark matter continues to col- 
lapse and virializes at the virial radius r^ii < rsh- Therefore, 
when we solve for hydrostatic equilibrium in the post-shocked 
gas, we assume that the dark matter has relaxed and we ex- 
trapolate the NFW profile between the virial- and shock radii. 
This inconsistency in our treatment of the gravitating mass is 
an unavoidable consequence of applying a quasi-equilibrium 
model for the dark matter outside the virial radius. 



4.1.3 Gas Equations 

We solve for the state of the gas using a method very similar 
to that described in section 3, but we now solve for the shock 
radius (,{m) self-consistently. We discretize the halo formation 
into the accretion of a finite number of shells and, for each 
shell rrii, we solve the full eigenvalue problem for its shock 
strength ^(mi). Thus, we simultaneously build up solutions for 
the shock radius and for the temperature profile as functions 
of time. As discussed above, we assume that the gas and dark 
matter first separate at the virial shock. The baryonic accre- 
tion rate at the shock radius is thus also proportional to 
equation 12. This is qualitatively consistent with the findings 
of Faucher-Giguere et al. (2011), who show that the baryonic 
accretion closely tracks the dark matter accretion history in 
simulations of high-mass halos. 

The equations for mass conservation and hydrostatic equi- 
librium are: 

;2 / r^./.,^\ 3/5 

(13a) 



dx 
dy 

dPi 
dy 



fKiiy) 



Pi{y). 

^2/3^8/3 I lQg(;^ ^ , 



c) — cx/{l + cx) 



X log(l -f c) — c/(l + c) 
with the boundary conditions: 



(13b) 



Pi{y = 


1) 


4^2 

^ ~3~ 


x{y = 


1) 


= 2^ 


x{y = 


0) 


= 0. 



2/3,8/3 



1-C 
(20^ 



(14a) 

(14b) 
(14c) 

In the above, rj = fdyn(c'lii A/vir/9t) measures the accretion 
rate of the halo, and h = H/Ho = [flm{l + zf + QaY^'^ is the 
Hubble parameter in units of Hq. In equations 13 and 14, 77, 
h, and ^ are all evaluated at the epoch of the most recently 
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accreted shell. It is straightforward to show from equation 12 
that m = my{l + z — b)H . Thus, the dimensionless accretion 
rate 77 is given hj rj = 7(1 + z — b)/10. 

The shock entropy Kf^ takes the form: 



y 



-1 2/3 



ay)] 



4/3 



(15) 



_v[4y)] X h[z{y)]_ 

where z{y) represents the redshift at which the baryonic shell 
y accreted and ^{y) represents the shock radius of the cluster 
at that redshift. 

Inserting the conductivity appropriate for a fully-ionized 
hydrogen plasma (Spitzer 1962) into equation 2 yields: 

47r ^^/^ ^ j^^ 



dlnKi 
dm 



0.70/sp 



vir, 



/b hmr] \10^^ Mq 



X 



5/2 



2 dTi 



+ 



2Ti 



dx 



+ 



9a; 2 



(16) 



(Recall that the coordinate m tracks the formation of the 
cluster and thus functions as a time coordinate). We integrate 
this equation between accretion events with an explicit, sub- 
cycled method. The conductive heat flux must vanish at the 
origin by spherical symmetry; thus, we adopt dT/dr = as 
the inner boundary condition for equation 16. The precise 
boundary condition on the heat flux at the shock is uncertain 
because it depends on the physics of coUisionless shocks in 
the presence of strong thermal conduction. We proceed by 
assuming that electrons do not diffuse across the shock into 
the upstream flow; in our model, the shock thus serves as an 
insulating boundary. This assumption is convenient because 
thermal conduction does not modify the structure of the shock. 
We have also tried calculations in which we keep the heat 
flux constant at the shock; the results with this alternative 
boundary condition were very similar to those we present here. 

Note that equations 13-15 are independent of halo mass, 
while equation 16 is not. Conduction thus introduces non-self- 
similar behavior and may influence mass-observable relations. 
We quantify this departure from self-similarity in the following 
section. 

The post-shock entropy (eq. 15) and the outer boundary 
condition on the pressure (eq. 14a) depend on the accretion 
history through the dimensionless accretion rate (cf. § 3). 
Thus, the diversity in accretion histories may generate scatter 
in the ICM properties at fixed halo mass. In order to estimate 
the statistics in ICM properties, we generate an ensemble 
of accretion histories for each halo mass Mvir, 0, with the 
exponents 7 and 6 drawn from the distribution in the appendix 
of McBride et al. (2009). This ensemble yields information 
about the statistics of the cluster population and the extent 
to which the variation in accretion histories creates scatter 
in the temperature profile and mass-observable relations. We 
present these results in the next section. 

4.2 Results 

Before studying the effect of thermal conduction, we present 
adiabatic models in which we set the effective conductivity 
to zero. These profiles facilitate comparison with the simpler 
models in section 3 and also provide a basis for understanding 
the models with conduction. 

Figure 5 shows representative, adiabatic solutions for the 



redshift-zero temperature profiles resulting from several differ- 
ent accretion histories. (Because we have neglected thermal 
conduction, these solutions are independent of halo mass.) In 
all cases, the temperature approaches the virial temperature 
near the scale radius of the halo, and falls by a factor of ~ 2 
by the virial radius. We note that the temperature profiles 
in figure 5 are not monotonic with density and thus cannot 
be described with polytropic models. Interestingly, however, 
a polytrope with an index 7 ~ 1.2 (as assumed by V03) pro- 
vides a reasonable fit between the scale radius and ~ 0.5rvir. 
These profiles are qualitatively very similar to the nfw models 
in figure 3 with tdyn/io ~ 0.15, as is reasonable for clusters 
(§ 3.2). The location of the accretion shock is also consistent 
with our interpretation in section 3: Type III halos, which 
experience slower accretion at late times, have larger shock 
radii than the more rapidly accreting Type halos. Thus, the 
intuition we developed in section 3 likely holds even for the 
more complex models in this section. 

Ignoring thermal conduction is not a well-motivated ap- 
proximation, however: as discussed in section 1, the timescale 
for heat to diffuse through massive clusters (r^i^/Xc ~ 1 Gyr) 
is shorter than the typical age of the ICM (~ 5 Gyr). Conse- 
quently, non-cosmological simulations of isolated halos (e. g. 
Parrish et al. 2008) show that the ICM becomes almost com- 
pletely isothermal after ~ 2 Gyr. By analogy, one might there- 
fore expect conduction to significantly modify the temperature 
profile shown in figures 3 and 5. 

Figure 6 compares the temperature profiles of clusters 
with different effective conductivities. In order to maximize 
the influence of thermal conduction, we show lO^^Af© clusters 
(which are hotter, and thus more conductive than lower mass 
clusters), with Type III accretion histories (which formed 
comparatively early and thus provide more time for conduction 
to operate). As expected, thermal conduction smooths out 
the temperature profile in the ICM; the effect, however, is 
substantially weaker than has been found in non-cosmological 
simulations (see above). The profiles we obtain do not become 
isothermal, in qualitative agreement with x-ray observations 
of the ICM (George et al. 2009; Simionescu et al. 2011), and 
also with cosmological simulations of clusters which include 
thermal conduction (Dolag et al. 2004; Burns et al. 2010). 
The disagreement with simulations of isolated halos suggests 
that cosmological accretion and the continued formation of 
clusters is essential not only for understanding the origin of 
the large-scale temperature gradient in the ICM (§ 3), but also 
for how this gradient persists in spite of thermal conduction. 

Two important effects differentiate cosmological calcula- 
tions of thermal conduction from non-cosmological simulations 
of isolated halos. First, even though the present-day timescale 
for thermal conduction is shorter than the mean age of massive 
clusters, the gas near the accretion shock is always young and 
thus unaffected by conduction. The jump conditions at the 
virial shock therefore still determine the outer temperature, 
even when we take thermal conduction into account. The sec- 
ond effect differentiating cosmological and non-cosmological 
calculations is that the ICM at higher redshift had a lower tem- 
perature and thus a much lower conductivity than it has at 
redshift zero. Even at small radii, the age of the ICM (~ 5 Gyr) 
thus over-estimates the timescale over which thermal diffu- 
sion operates. These two effects, which are not present in 
non-cosmological simulations, strongly limit the influence of 
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Figure 5. Representative temperature profiles at redsliift z = 
generated by eacli type of accretion iiistory from McBride et al. 
(2009) (see fig. 4). The profiles are calculated using the Lagrangian 
method from section 4, but with no thermal conduction (the 
temperature profiles are thus independent of the halo mass). Note 
the qualitative agreement with the simple models from section 3. In 
addition, the location of the shock radius depends fairly sensitively 
on the accretion history of the halo. 
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Figure 6. Influence of the effective conductivity on cluster temper- 
ature profiles. These curves show temperature profiles for massive, 
10^^ Mq halos with Type III accretion histories at redshift 2 = 0. 
This choice of mass, accretion type, and redshift maximizes the 
effect of thermal conduction; nonetheless, the change to the tem- 
perature profile at large radii is modest. As discussed in the text, 
the marked effect of thermal conduction on the temperatures at 
small radii is not a firm prediction, since our models neglect both 
radiative cooling and feedback heating, which are important in 
this region (Voit 2011). 



thermal conduction on the large-scale temperature profile of 
the ICM. 

We illustrate these points in figure 7, where for each 
model from figure 6 we plot the ratio of the age of each shell 
in the ICM to its field-free conduction timescale. We define the 
conduction timescale as icond = \d\nK/dt\~^ , calculated using 
equation 2 with fsp = 1. We define the age of the shell tago 
as the time over which conduction operates at its present-day 
efficiency; we approximate this as the minimum of the time 
since the shell accreted and \d\mn/dt\~^ , the timescale over 
which the ICM conductivity changes appreciably. Conduction 
should limit the temperature profile of the ICM to everywhere 
satisfy fagc/icond ^ fsp- Thus, by looking at the adiabatic 
models in figure 7 (blue, with fsp = 0), we can estimate how 
large of an effect conduction would have if we were to include 
it. 

Figure 7 shows that the ratio tago/icond falls steeply near 
the virial shock, again reflecting that the recently accreted 
gas is too young to be influenced strongly by conduction. At 
smaller radii, however, tage/icond ~ 5 in the adiabatic model, 
suggesting that thermal conduction plays a significant role in 
determining the temperature profile in the ICM. The thin blue 
curves in figure 7 show that this conclusion depends strongly 
on redshift, however. The flattening of the temperature proflle 
shown in flgure 6 is thus a fairly recent phenomenon and 
is far less pronounced at redshift z > 0.2. Recall also that 
flgures 6 and 7 show results for massive, 10^^ clusters with 
Type III accretion histories. Clusters with lower masses or 



different accretion histories are likely to be even less strongly 
influenced by conduction (see flgure 9, below). 

The red and yellow curves in figure 7 show that when we 
include thermal conduction with a given suppression factor 
/sp, the temperature gradient adjusts so as to keep the ratio 
icond/iagc below fg^- The ratio tcond/iagc is proportional to 
the gradient term; 

"2 dTi J_fdT\y d'^Ti ' 
x~dx ^ 2I\ \~dx ) dx^ ■ 

Thus, thermal conduction can either make the ICM isothermal 
(so that the heat flux vanishes), or make T ~ r"^''^ (so that 
the heat flux is constant). We flnd that our model clusters 
initially take the second approach; the temperature proflle 
adjusts so as to have a nearly constant heat flux as a function 
of radius. This is not possible near the origin, however, where 
spherical symmetry requires the heat flux to vanish. Thus, 
over a longer timescale, the center of the ICM cools and the 
ICM begins to become isothermal. The yellow curve in figure 6 
shows the beginning of this process, but there is not enough 
time before redshift z = Q for a signiflcant fraction of the 
ICM to become completely isothermal, even in cases where 
conduction is most effective. 

Though conduction may have considerable effect on the 
temperature proflles of massive clusters by redshift zero, its 
influence on their entropy profiles is less pronounced. We show 
this in flgure 8, which compares entropy proflles for massive, 
clusters with and without conduction. Between the 
scale radius and the virial radius, both models agree fairly 
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Figure 7. Ratio of the age of each shell in the ICM to its field- 
free conduction timescale (see text for definitions). The models are 
the same as in figure 6. This ratio is defined using the field- free 
conductivity (i.e. for fsp = 1). Models with finite suppression 
factors should satisfy tage/tcond $ fsp ' thick gray lines show 
this limit for /gp = 1 and 1/3. Adiabatic models predict tcond $ 
(0.2 — 0.3)tage at intermediate radii; as a result, conduction alters 
the temperature profile and maintains <agc/tcond J: fs„ ■ 




10-1 1 

Figure 8. Entropy profiles in calculations with different values of 
the conductivity. As in figure 6, this figure shows 10^^ Mq halos 
with Type III accretion histories; these properties maximize the 
effect of thermal conduction. Also shown is the K ^ r^'^ power-law 
derived in V03. 



well with a power-law fit derived from x-ray observations 
(e.g. Cavagnolo et al. 2009). Interestingly, tliermal conduction 
smooths out some of tiie non-power-law behavior introduced 
by the accretion history, leading to an entropy profile more 
similar to the adiabatic self-similar profile. 

Figures 5-8 illustrate the effects of cosmological accretion 
and thermal conduction on the temperature of the ICM in 
individual clusters. As mentioned in section 4.1, in order to 
study the statistics of temperature profiles as a function of 
virial mass, we run an ensemble of models with accretion his- 
tories drawn from the probability distribution in the appendix 
of McBride et al. (2009). The upper panels of figure 9 show 
these temperature profiles for three different halo masses with 
(yellow) and without (blue) thermal conduction. In the models 
with conduction, we have assumed a conductive suppression 
factor /sp = 0.3. The temperature profile roughly follows the 
pattern found in section 3, but with fairly significant, ~10%, 
scatter. This scatter is caused purely by the variation in accre- 
tion histories. We have neglected several processes, including 
mergers, heating by dynamical friction, radiative cooling, and 
heating by AGN outflows, all of which may increase the scatter 
above that shown here. 

Figure 9 also illustrates the mass-dependence of conduc- 
tion's role in the ICM. Because the thermal conductivity of a 
plasma depends sensitively on its temperature as k oc T^/^^ 
higher-mass halos are more strongly influenced by conduction. 
More quantitatively, the conduction timescale in the ICM falls 
with mass: tcond ~ ''vir/Xc ~ l/A/vir- Figure 9 shows that the 
normalization is such that conduction has a minor influence 
on W^^Mq halos, ^ but that it is significant for more massive, 
IO^^Mq clusters. 

Finally, in the upper-right panel of figure 9, we have over- 
laid observational data from Leccardi & Molendi (2008) and 
Simionescu et al. (2011).^ The data from Leccardi & Molendi 
(2008) (shown in black) represent an average temperature pro- 
file derived from 48 clusters from the XMM-Newton archive. 
The points from Simionescu et al. (2011), shown in green, 
represent suzAKU observations of the Perseus cluster out to 
r2oo- Outside the scale radius, where our method is applicable, 
our model agrees favorably with the data. Though this agree- 
ment is encouraging, we caution that both the normalization 
and slope of the data points in this plot are sensitive to the 
assumed virial mass of the halo, which is uncertain in x-ray 
observations. Furthermore, several uncertainties in our model, 
including the turnaround radius, the effect of non-smooth 
accretion, and deviations from spherical symmetry preclude a 
very quantitative comparison with the data. 

4.3 Mixing-Length Theory and the MTI 

As an example application of our models, we use mixing-length 
theory to estimate the turbulent pressure support produced 
by convection in clusters. In dilute, magnetized plasmas such 
as the ICM in galaxy clusters, convective stability depends 

This conclusion applies to the large-scale temperature gradient in 
the ICM. Of course, small-scale features may be strongly infiuenced 
by conduction, even in IO^'^Mq halos (cf. Dolag et al. 2004). 
^ The observation presented in George et al. (2009) provides another 
useful constraint on the properties of the ICM at large radii. We do 
not include it in figure 9, however, because it is not de-projected 
and thus not directly comparable to our results. 
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Figure 9. (Top row): Temperature profiles with (yellow; /gp = 0.3) and without (blue) thermal conduction for different halo masses. Black 
lines show the median profile from our ensemble of accretion histories based on figure 1 and colors show contours enclosing 50%, 90%, and 
99% of the models. The points with error bars in the top-right panel show temperature profiles derived from x-ray observations. The black 
and green points are taken from Leccardi &c Molendi (2008) and Simionescu et al. (2011), respectively. (Bottom row): Estimated fraction of 
non-thermal pressure support generated by the MTI. See section 4.3 for details. 



on the temperature gradient of the plasma (Balbus 2001, 
later generalized by Quataert 2008 and Kunz 2011). This 
convection, known as the mti, may produce strong turbulence 
in clusters (McCourt et aL 2011; Parrish et al. 2012) and 
thus may provide enough non-thermal pressure support to 
bias hydrostatic mass estimates of cluster halos. Parrish et al. 
(2012) found that the convective velocities produced by the 
MTI roughly obey mixing-length theory. This motivates us to 
use our model temperature profiles to estimate the turbulent 
pressure support produced by the mti as a function of cluster 
mass and redshift. We note, however, that Kunz et al. (2012) 
have shown that the strength of this turbulence depends on 
the magnetic field strength in the ICM, an effect which we do 
not account for in our simple estimates. 

Assuming that the convective motions retain their co- 
herence for a fraction a of a pressure scale-height and that 
magnetic tension does not suppress the convective motions, 
we expect the instability to drive turbulent convection with 
Mach numbers of order: 

M^a^H — j , (17) 



where H = {d\n P/dr)~^ is the pressure scale-height. The 
bottom panels of figure 9 show this estimate for the turbu- 
lent pressure support (proportional to M^). Assuming (as 
suggested by the sinralations in Parrish et al. 2012) that the 
mixing length parameter a ~ 0.5, this figure suggests that 
turbulence driven by the mti contributes ~ 5 — 10% of the 
pressure support outside of rsoo- 

Since the mti is not suppressed by other sources of turbu- 
lence (McCourt et al. 2011; Parrish et al. 2012), this pressure 
support adds to that already present due to turbulence driven 
by infalling substructure or by galaxy wakes. This analysis 
suggests that the mti plays an important role in the dynamics 
of the ICM, especially at radii > O.Srvir. Thus, the mti may 
provide an interesting correction for hydrostatic mass esti- 
mates of cluster halos. Unfortunately, both the large scatter 
in the strength of the mti and its strong radial dependence (cf. 
Parrish et al. 2012) seem to preclude a simple fitting function 
for the fraction of turbulent pressure support as a function of 
halo mass. 

By applying mixing-length theory in equation 17, we 
have implicitly assumed that the mti grows rapidly enough 
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to establish convection by redshift zero. We find that tuTi ^ 
0. Stage within rvir; thus, the Mach numbers in figure 9 (and 
also in Parrish et al. 2012) are likely to be reasonable estimates. 

Another assumption implicit in our use of mixing-length 
theory is that conduction is rapid enough to sustain the mti. 
Thus, the results in the lower panels of figure 9 are only 
valid when the conduction timescale across an unstable mode 
~ (aH)^ /x is less than the growth time of the mti. Massive 
clusters (Mvir lO^**'^ Mq), likely satisfy this ordering of 
timescales, but lower masses halos (Mvir < 10" Mq) may 
not. Thus, the mass-dependence of the mti is not likely to 
be monotonic. In the most massive halos (Mvir lO'^^ Mq), 
thermal conduction is more efficient and weakens the temper- 
ature gradient at intermediate radii (fig. 9). In lower mass 
halos (Mvir 10^'' Mq), on the other hand, conduction may 
not be fast enough to drive the mti to its full potential. This 
analysis suggests that lO^**'^ Mq halos experience the most 
vigorous convection driven by the mti. 



5 DISCUSSION 

This paper provides a simplified, spherically symmetric model 
for the temperature profiles of the hot plasma in galaxy groups 
and clusters. Our model is similar in spirit to earlier studies of 
the entropy profiles in clusters (e. g. Tozzi & Norman 2001; Voit 
et al. 2003), but builds on these earlier studies by focusing on 
temperature and by including the effects of thermal conduction. 
Our results agree reasonably well with the profiles derived 
from x-ray observations (cf. Leccardi & Molendi 2008) and 
from numerical simulations (e.g. Dolag et al. 2004). 

We have shown that the large-scale temperature gradient 
in the ICM is primarily determined by the accretion history 
of its halo: while the gas near the center of a cluster reaches 
the virial temperature of the halo, the temperature at the 
virial shock is determined by the ram pressure of the accreting 
gas. This difference sets the overall shape of the temperature 
profile. The timescale for thermal conduction (tcond ~ f'vir/Xo) 
is somewhat shorter than the age of the ICM. However, the 
influence of thermal conduction on the global temperature 
profile in clusters is mitigated for two reasons. The gas near 
the center of the cluster is less strongly effected by conduction 
because it was cooler in the relatively recent past. The gas 
near the virial shock, on the other hand, has only recently 
accreted and is younger than the conduction timescale. Thus, 
conduction has a diminished effect both near the center of 
the cluster and near the outskirts. As a result, it does not 
dramatically change the mean temperature profile (fig. 9). 
Of course, conduction can have a dramatic effect on small- 
scale inhomogeneities in the ICM (cf. Dolag et al. 2004); such 
inhomogeneities cannot be studied in our one-dimensional 
model. 

Our results demonstrate the close relationship between 
the temperature gradient in clusters and the cosmological 
evolution of the host halo. This implies that numerical studies 
of isolated cluster models (e. g. Parrish et al. 2008) cannot 
correctly predict the evolution of the large-scale temperature 
profile, though they are very useful for studying other aspects 
of the ICM, such as the interplay among cooling, feedback, and 
plasma instabilities within the scale radius of the halo. 

One of our motivations for studying the effects of conduc- 
tion and halo accretion history on temperature gradients in 



clusters is that the free energy in the non-zero temperature 
gradient drives an efficient convective instability, the mti. The 
results in section 4.3 (e. g. fig. 9) show that the turbulent pres- 
sure support generated by the mti may be of order ~ 5 percent 
of the thermal pressure, and that it scales non-monotonically 
with halo mass. The magnitude of the turbulent pressure sup- 
port is sensitive to the accretion history and does not seem 
amenable to a simple fitting formula. 

Halo accretion histories have been studied extensively 
with numerical simulations. We use the fits to the Millennium 
simulation from McBride et al. (2009) to estimate the scatter 
in temperature profiles as a function of halo mass; at redshift 
z — 0, this scatter is of order 10%. This scatter likely con- 
tributes to the dispersion in cluster mass-observable relations 
relevant to x-ray and SZ observations. Perhaps more interest- 
ing are the effects of thermal conduction and convection, which 
introduce systematic changes to the temperature profile with 
mass. In particular, conduction smooths out the temperature 
profile (and decreases the peak temperature in the halo) by an 
amount that increases monotonically with halo mass (fig. 9); 
convection, on the other hand, produces turbulent pressure 
support that is non-monotonic in halo mass, peaking around 
lO^-^Afe halos (§ 4.3). 

Figures 5, 6, and 9 demonstrate that the effect of ther- 
mal conduction on a cluster's temperature profile is at least 
as large as the differences produced by normal variation in 
accretion histories. Any variation in the factor fsp (which 
parameterizes the suppression of the effective radial thermal 
conductivity relative to the field- free value), if it exists, would 
create additional scatter in the temperature profiles at fixed 
mass. Possible effects infiuencing fsp include magnetic draping 
around infalling or orbiting substructure (Dursi & Pfrommer 
2008; Pfrommer & Dursi 2010) and the strength of the mag- 
netic field. These processes may also contribute to the scatter 
in cluster mass-observable relations. 

The models presented in this paper provide a simple 
explanation for the physics that sets the temperature profiles 
in galaxy groups and clusters at large radii. Our results are 
consistent with current observational constraints on cluster 
temperature profiles at large radii (fig. 9). They also highlight 
several processes which may bias hydrostatic mass estimates 
of clusters, including modifications to the scaling relation 
T(M) by conduction and by turbulence driven by the mti. 
The approximations in this paper, especially the assumptions 
of smooth accretion and spherical symmetry, preclude precise 
estimates of the non-thermal pressure support produced by 
the MTI. These limitations can be addressed using cosmological 
simulations. 
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